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Abstract 

The continuum shape of 26 broad absorption hne (BAL) QSOs is shown to be virtually identical to 
that of non-BAL QSOs when the spectra of the former are dereddened. A procedure for dereddening of 
BAL QSO spectra is introduced, based on measurement of the flux in three bands and determination of 
the colour excess in comparison with an unreddened spectrum. This method allows us to derive the type 
and amount of extinction. We show that reddening in high-ionization BAL QSOs is comparatively small 
and the range of extinction laws varies from Milky Way, LMC to SMC type. On the other hand, low- 
ionization BAL QSOs seem to experience reddening which can be attributed to SMC type extinction alone. 
Comparison of the colour indices of composite spectra from Francis et al. (1991) and Zheng et al. (1997) 
with those derived from individual BAL QSO spectra shows that the Francis et al. (1991) composite might 
be a more suitable representative of an unreddened spectrum. In general, we conclude that the spectral 
energy distribution of normal QSOs can be used as an incident continuum shape for BAL QSOs. 

Key words: Interstellar: dust — Interstellar: extinction — Quasars: general — Quasars: absorption 
lines — Ultraviolet: galaxies 
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1. Introduction 

In optically selected QSO samples, ~10% are broad ab- 
sorption line (BAL) QSOs (Weymann et al. 1991). BAL 
QSOs can be subdivided into high- and low-ionization 
(hi-BALs and lo-BALs) where hi-BALs have absorptions 
from high ionization species (Nv, Civ, and Siiv) and lo- 
BALs, which make up only ~10% of the BAL QSO pop- 
ulation, have additional absorptions from low-ionization 
species (Aim, Alii, and sometimes Mgii). In some cases 
lo-BALs also have Fell and Fein absorptions (iron lo- 
BALs). 

Becker et al. (1997) reported the discovery of two new 
iron lo-BALs from radio-selected samples (FIRST survey; 
Becker et al. 1995) with one of the objects (1556-1-3517) 
showing severe extinction in the spectrum. Goodrich 
(1997) pointed out that a significant number of BAL 
QSOs are missed in optically selected samples due to at- 
tenuation of the continuum and that the true fraction of 
BAL QSOs might be as high as 30% instead of 10%. This 
could mean that these heavily reddened objects are more 
common among QSOs than previously thought. 

The continuum shape of BAL QSOs is often flat or dis- 
torted compared to that of non-BAL QSOs (Weymann 
et al. 1991). Sprayberry and Foltz (1992) showed that 
the composite spectra of high- and low-ionization BAL 



QSOs can be brought into coincidence by dereddening the 
lo-BAL composite spectrum relative to the hi-BAL com- 
posite with appropriate extinction laws. This method 
shows only the tendency for lo-BALs to be redder than 
hi-BALs on average. 

The present work describes a systematic method for es- 
timating the amount of extinction for individual hi- and 
lo-BAL spectra. The systematic method for deriving ex- 
tinction for all BAL QSO types allows a homogeneous 
comparison of the whole data set, which can also be ex- 
tended to non-BAL QSO spectra. Section 2 describes 
the observational data and section 3 will introduce the 
dereddening procedure for BAL QSOs followed by results 
and discussions in section 4. Conclusions arc shown in 
section 5. 

2. Observational Data 

The observational data consist of a total of 26 BAL 
QSO spectra; 15 objects are from Weymann et al. (1991), 
8 from Korista et al. (1993), 2 from Brothcrton et al. 
(1997), 1 from Ogle (1997, private communication) and 
parts of spectra from Becker et al. (1997) and Wills (1998, 
private communication) have kindly been made available 
by the authors. From the sample we chose only those 
spectra which had spectral coverage of all three bands 
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Table 1. Main parameters of the program quasars. 



Object 



BAL type 



F1750/F: 



1750/ -r 2 100 



F2100/F:-. 



2100/ -f^ 3000 



Extinction type 



Ref. 



0004+0147 

0018+0047 

0025-0151 

0051-0019 

0054+0200 

0059-2735 

0145+0416 

0840+3633* 

0932+5010 

1009+0222 

1205+1436 

1208+1535 

1212+1445 

1216+1103 

1228+1216* 

1239+0955 

1303+3080* 

1331-0108* 

1556+3517* 

1641+4115 

2154-2005 

2201-1834 

2202-2007 

2225-0534* 

2350-0045 

2358+0216 



1.706 


hi-lo 


1.816 


hi 


2.072 


hi 


1.705 


hi 


1.868 


hi 


1.584 


iron 


2.029 


hi 


1.225 


iron 


1.914 


hi-lo 


1.343 


hi-lo 


1.629 


hi 


1.956 


hi 


1.621 


hi 


1.615 


hi 


1.408 


hi 


2.014 


hi 


1.760 


hi 


1.867 


lo 


1.487 


iron 


2.005 


hi 


2.028 


hi 


1.817 


lo 


2.188 


lo 


1.981 


hi-lo 


1.626 


hi-lo 


1.854 


lo 



1.16 
1.37 
1.25 
1.23 
1.28 
1.09 
1.22 
1.01 
1.30 
1.11 
1.27 
1.37 
1.22 
1.41 
1.32 
1.27 
1.42 
0.88 
0.60 
1.25 
1.38 
0.94 
1.21 
1.15 
1.28 
0.78 



1.68 
1.72 
1.51 
1.31 
1.39 
1.33 
1.44 
1.16 
1.75 
1.57 
1.46 
1.44 
1.44 
1.42 
1.44 
1.96 
1.47 
1.38 
0.27 
1.45 
1.68 
1.10 
1.83 
1.09 
1.53 
1.08 



sihcate/SMC 


b 
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SMC-LMC 
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SMC-LMC 
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LMC 
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SMC 


a,e 


SMC-LMC 
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SMC 


c,d,f 
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sihcate/SMC 
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SMC-LMC 


a 


LMC-MW 
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SMC-LMC 
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LMC-MW 
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LMC 
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LMC-MW 
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sihcate 
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SMC 


c 


SMC-LMC 


a 


- 


a 


sihcate 


a 


- 


b 


SMC-LMC 


a 


SMC-LMC 


a 


silicate 


b 



a: Weymann et al. (1991); b: Korista et al. (1993); c: Brotherton et al. (1997); d: Becker et al. (1997) e: Ogle (1997, 
private communication); f: Wills (1998, private communication). 
* Spectrum shown in figure 5. 



defined below in order to do the colour measurements. 

We divide the quasars into hi- and lo-BAL types (hi 
and lo in table 1), similar to the classification in Wey- 
mann et al. (1991). However, the classification for hi- 
and lo-BALs is somewhat ambiguous and for those cases 
where it is not clear which group an object belongs to 
we indicated this with "hi-lo" . Additionally we use the 
classification of the iron lo-BALs (iron in table 1) which 
was introduced by Becker et al. (1997). Table 1 (columns 
1, 2, and 3) lists the names of the program BAL QSOs, 
the redshifts and the classifications of the objects. 



3. Dereddening Procedure 

The extinction laws used here are the Milky Way 
(MW), Large Magellanic Cloud (LMC), and Small Mag- 
ellanic Cloud (SMC) from Pei (1992), the attenuation 
curve from Calzetti et al. (1994) and theoretical extinc- 
tion curves for astronomical silicate with different grain 
sizes from Laor and Draine (1993). 

3.1. The Unreddened Continuum 

Sprayberry and Foltz (1992) established the idea for 
the existence of dust extinction in BAL QSOs by show- 



No. 3] 



On the Continuum Shape of BAL Quasars 



ing that lo-BALs are reddened relative to hi-BALs. We 
take a more general approach by comparing BAL QSOs 
with a representative non-BAL QSO sample. By com- 
paring BAL QSO spectra with a single spectral shape we 
assume that all BAL QSOs are reddened relative to the 
representative spectrum, which can be understood as an 
"intrinsic quasar spectrum". We refer to subsection 4.1 
for discussion of the intrinsic quasar spectrum. 

For the present investigation the following spectra are 
taken into consideration for determining the continuum 
shape of an unreddened QSO spectrum: (a) the compos- 
ite spectrum from Francis et al. (1991; hereafter FSED), 
(b) the composite spectrum Zheng et al. (1997; hereafter 
ZSED), and (c) an artificially steepened FSED (hereafter 
FSED-plus) with the spectral index (defined by F^, ~ j/") 
corresponding to a = — 1 between 1250 and 2100 A and 
a = -1-0.4 between 2100 and 3100 A. We note that spec- 
tra of BAL QSOs are included in the composite spectrum 
from Francis et al. (1991) but the final spectral slope 
of the composite was chosen to have the median slope 
and thus the BAL continuum will hardly affect the final 
shape. 

The spectral slopes of FSED-plus result from the fol- 
lowing consideration: A comparison (e.g. figure 9 in 
Zheng et al. 1997) shows that FSED and ZSED have an 
almost identical continuum shape between 1250 A and 
2100 A with a « —1. However, longward of 2100 A the 
Francis et al. (1991) composite spectrum is harder with 
a K, —0.3, whereas the Zheng et al. (1997) remains flat 
(a K, —1). Examination of the distribution of spectral 
slope (Francis et al. 1991, figure 1) shows that a has a 
wide spread around the median value of a = —0.32. Thus 
we can consider a = —1 (soft spectrum) as a lower bound 
and a ~ +0.4 (hard spectrum) as the upper bound for 
the spectral index for wavelengths > 2100 A. Our choice 
of FSED-plus together with the composite spectra FSED 
and ZSED should cover most of the range in the spectral 
index. 

For the sake of clarity later in figure 5 we approximate 
the continuum shape of the composite spectrum FSED 
as a broken power-law with a = — 1 for 1200 A < A < 
2100 A and a = -0.32 for A > 2100 A. Note that all 
measurements are carried out on the composite spectra 
and not on the approximate power-law. 

3.2. Definition of the Bands 

We define three bands in the spectrum at restframe of 
the quasar around 1750 A, 2100 A, and 3000 A. Mea- 
surement of the continuum level of BAL QSO spectra 
is difficult, since BAL spectra have numerous absorption 
lines of high- and low- ionization elements that arise in the 
broad and narrow line region of the quasar. Therefore, 
windows, which seem relatively free from either emission 
or absorption, were defined within the bands as follows: 
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Fig. 1. The windows in the bands B1750, B2100, and B3000 are 
shown. The two spectra shown are 0840+3633 and 2225—0534. 



B1750: at 1715-1735 A and 1760-1780 A 

B2100: at 1980-2010 A, 2125-2145 A, and 2185-2195 A 

B3000: at 2880-2930 A and 3030-3090 A 

The total coverage of the bands is AA w 1250 A with 
separation of 350 A and 900 A between the bands. This 
is sufficient to trace extinction effects in the UV region 
where the extinction curves rise steeply. However, the 
necessity for a wide spectral coverage restricts the num- 
ber of objects which can be analyzed. Note that B2100 
is located near to the 2175 A feature of the Galactic ex- 
tinction curve. 

Figure 1 shows the position of the windows in the in- 
dividual bands. As an example we show two BAL spec- 
tra, one which has relatively strong absorption lines of 
low and high ionization species (lo-BAL 0840+3633) and 
another, which has only high ionization lines (hi-BAL 
2225-0534). 

The choice of the windows can be summarized as fol- 
lows: 

Band around 1750 A: (figure 1, upper panel) This 
band lies between C iv and the semiforbidden emission 
line of Cm]. Two windows (la) at 1715-1735 A and 
(lb) 1760-1780 A cover a total of 40 A. The strong ab- 
sorption lines are the singlet AliiA1671, the doublets 
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Fig. 2. The effect of extinction on the colour indices is shown, (a) for the case of the artificially steepened Francis ct al. (1991) 
spectrum FSED-plus defined in subsection 3.1 (b) for Francis et al. (1991) composite spectrum and (c) for Zheng et al. (1997) 
composite spectrum. The diff'erent lines show the effect of extinction when the empirical extinction laws of the MW, LMC, and 
SMC arc applied. The type of extinction law is indicated in the figure. The extinction increases along the lines starting from the 
filled circle at the top of the lines which represents the measurement of the unreddened composite spectrum. The observations (open 
circles) are shown with errorbars. The dots on the lines mark the amount of extinction increasing by a factor of 2 from dot to dot. 
The thick dot on each line indicates E{B — V) = 0.25. 



Si II AA 1808, 1817, and Al ill AA 1855, 1863. The Si 11 line 
is only AA = 28 A away (which corresponds to a veloc- 
ity of w = cAA/A ~ 4 000 kms"^) from window (lb) 
but the blueshift of broad absorption lines of Sin sel- 
dom exceed velocities of 4000 kms^^ and thus are un- 
likely to affect the measurement in this window. Al ill 
shows much stronger BALs. However it is unlikely for 
the width of the absorption line to exceed 13 000 kms^^ 
which would diminish measurements in (lb). The weaker 
features are absorptions from multiplet lines of Fe 11 and 
possibly Ni ill. In strongly absorbed quasar spectra these 
weak emissions and absorptions may have an effect on 
window (la). The emission-like feature around 1790 A is 
caused by the Fell uv 191 transition. 

Band around 2100 A: (figure 1, middle panel) The 
second band lies between Cm] and Mgii. Three win- 
dows at 1980-2010 A, 2125-2145 A, and 2185-2195 A 
cover a total of 60 A. This definition is close to the one 
used by Weymann et al. (1991). Window (2a) is posi- 
tioned before the rise of the semiforbidden line emission 
of C III] A 1909 on the red side. The soft rise between (2a) 
and (2b) with its peak around 2070 A is likely caused by 
emission from the Fe II multiplet. This iron multiplet also 
shows BALs in some quasars but its absorption is unlikely 
to reach window (2a) which is 9 000 kms~^ away. The 
well-known strong multiplets of Fen uv 1, 64 (around 
2600 A) and Fell uv 2, 3, 35, 36 (around 2400 A) pro- 
duce deep blueshifted BALs in 0840-f 3633. The series of 
three strong absorptions is likely caused by Fe n (around 
2250 A) and possibly Niii (around 2210 A and 2170 A). 
Windows (2b) and (2c) are set to avoid these spectral 
features. 



Band around 3000 A: (figure 1, lower panel) The third 
band is on the red side of Mgii. Two windows at 2880- 
2930 A and 3030-3090 A cover a total of 110 A. The main 
feature here is the emission lines of MgiiAA 2796, 2803. 
The continuum level around Mgii is raised by numer- 
ous Fe II multiplets and the Balnier continuum of hydro- 
gen. Although the continuum to the red of Mg II seems 
smooth, there appears a weak emission feature around 
3000 A, possibly caused by Fen, which varies in strength 
from object to object. Window (3a) is placed between 
Mgii and this emission feature and (3b) is placed to the 
red side of the emission. 

3. 3. Calibration of the Diagnostic Diagram 

The measurement of the fluxes i^i75o, .f2i00j a-nd i^sooo 
in the bands is obtained by integrating the spectrum over 
the windows and dividing it by the total width of the 
bands, i.e. average flux per A. Although the windows 
have been carefully chosen, it is important to check e.g. 
for overlapping BALs when measuring the bands and to 
correct the measure if necessary. 

Table 1 (columns 4 and 5) gives the derived colour in- 
dices which are defined as i^i75o/.F2ioo and i^2ioo/-F3ooo- 
The two-colour diagnostic diagram (TCDD) is obtained 
by plotting colour indices as (x, y)-pairs. The errors 
(shown in figure 2) are estimated by dividing the maxi- 
mum difference of fiux in a band with the square root of 
the number of bins within the integration interval. 

The effects of reddening with the various extinction 
laws, including SMC, LMC, and MW type from Pei 
(1992) and astronomical silicate from Laor and Drain 
(1993), are investigated by applying them to the com- 
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Fig. 3. Same as figure 2 but for astronomical silicate with dilTcrcnt grain sizes. The thick dots indicate E(B — V) = 0.1. Only the 
lines corresponding to the smallest and the biggest grain sizes are marked. For clarity, crrorbars have been omitted in this figure. 



posite spectra from Francis et aL (1991) and Zheng 
et aL (1997). After reddening the composite spectra the 
intensities in the bands are determined with the same 
procedure as for the observed spectrum. 

For simphcity, we use single grain size for the extinc- 
tion curve of astronomical sihcate where the foUowing 
sizes are used (given is the logarithm of grain sizes in 
Aim): -3, -2, -1.5, -1.35, -1.25, and -.75. 

Numerous niultiplet lines of iron-group elements can 
either depress the continuum or produce emission fea- 
tures in the bands and influence the measurements of the 
colour indices. Unfortunately, at this stage model calcu- 
lations of the low-ionization region cannot provide hints 
of whether the continuum level in the bands is lowered 
or raised. At present, our definition of the colour indices 
provides a homogeneous measurement for all quasars, in- 
cluding lo-BALs, and yields a robust estimate of the con- 
tinuum shape. 

4. Results and Discussions 

The extinction of the quasars can be obtained from the 
TCDD once the comparison spectrum for the diagram is 
chosen (see subsection 4.1). In order to distinguish the 
type of extinction more accurately we introduce a clas- 
sification scheme for the determination of the extinction 
type in subsection 4.2. The information obtained from 
the TCDD is used to dereddcn the spectra and we show 
examples for six of the program quasars at the end of 
this section. 

4-.1. The Two-Colour Diagnostic Diagrams 

Figure 2 shows the positions of the program quasars in 
the TCDD together with lines produced from FSED-plus 
(panel a), FSED (panel b), and ZSED (panel c) applying 
empirical extinction laws. 



Figure 3 shows the same program quasars as in figure 2 
but this time lines are produced by applying theoretical 
extinction curves of astronomical silicate with different 
grain sizes to FSED-plus, FSED, and ZSED (panels a, 
b, and c respectively). By comparing figures 2 and 3 
it becomes apparent that intermediate and large grain 
sizes mimic the lines in figure 2 that are produced by 
SMC and LMC extinction law respectively. The region 
above the line produced by SMC in panel b of figure 2 
may be explained by reddening FSED using extinction 
curves of astronomical silicate with small grain sizes. 

The first thing to notice in figure 2 is that the bulk of 
observations have i^2ioo/-p3000 values which are higher 
than those derived from ZSED at the top where the 
curves join (figure 2 panel c). We attribute this to a 
poor representation of the flux around 3000 A where the 
composite is derived from only a few quasars. The posi- 
tion of FSED in the TCDD, on the other hand, is at one 
of the bluest positions in the diagram (panel b). FSED- 
plus is further to the blue in the diagram (panel a) as 
expected since we chose a value resulting in the hardest 
(bluest) spectrum in our comparison. 

Nonetheless, the upper sequence of observations (see 
e.g. figure 2 panel b, to the left of the SMC-line) imphes 
the existence of a yet bluer quasar spectrum. Such an 
intrinsic spectrum would have a spectral slope which is 
harder than the Francis et al. (1991) composite. The arti- 
ficially hardened FSED-plus is only one way to represent 
the intrinsic spectrum, however, it should be reminded 
that FSED-plus is only a hypothetical spectrum whereas 
FSED and ZSED arc derived from observations and is 
therefore less desirable. Another method of obtaining 
the intrinsic spectrum of quasars is to take the steepest 
spectra from a quasar survey and construct a composite 
spectrum. If such an intrinsic spectrum exists it would 
mean that almost all quasar spectra observed are red- 
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Fig. 4. The effect of extinction on the Francis et al. (1991) composite spectrum for the empirical extinction curves MW, LMC, SMC, 
and attenuation curve (AC) (Calzetti et al. 1994) is shown. For the discussion is shown the line produced by astronomical silicate with 
grain size a = 0.032/im. The different symbols indicate hi-BALs (plus signs), lo-BALs (open circles), iron lo-BALs (open squares), 
and the quasars labeled as "hi-lo" in table 1 (open triangles). For comparison we have included the position of the average spectra 
from Weymann et al. (1991) non-BAL (filled diamond), lo-BAL (filled triangle), and hi-BAL (filled square). The filled circles on the 
lines mark the amount of extinction as described in figure 2. 



dened by dust. 

4.2. Determination of the Extinction Types 

From the discussion above we conclude that the ZSED 
and FSED-plus are less suitable for representing the in- 
trinsic shape of BAL QSOs. Therefore, the determina- 
tion of the amount and type of extinction will be carried 
out using FSED. 

Figure 4 shows the TCDD with lines for the case of 
Francis et al. (1991) composite spectrum when empiri- 
cal extinction laws are applied. From the position of 
the curves that are produced using the SMC, LMC, and 
MW extinction laws we define the following classification 
for the extinction type: SMC, SMC-LMC, LMC, and 
LMC-MW. Additionally we define the "silicate" type for 
objects on the left of the SMC-line. We subdivide the 
regions between SMC, LMC, and MW and decide upon 
the positions of the objects in the TCDD between the 
extinction types. In cases where it does not seem clear 
whether the spectrum has silicate or SMC type extinction 



this was marked with "silicate/SMC" in table 1. Note 
that i^2ioo/^300o niainly provides information about the 
amount of dust while -F1750/-F2100 generally distinguishes 
between extinction types. 

The lo-BALs with iron absorptions (0059—2735 and 
0840+3633; open squares in figure 4) lie on the SMC- 
line. Within its errors, even the iron lo-BAL 1556+3517, 
which suffers severe extinction, lies near the SMC-line 
(see panel b of figure 2). An alternative solution for 
iron lo-BALs gives astronomical silicate with grain size of 
about 0.05 yum. However, better agreement of the contin- 
uum level at shorter wavelengths between Civ A 1549 and 
Si IV A 1394 is achieved when the spectra are dereddened 
with the SMC extinction law. 

The lo-BALs without iron absorption (1331—0108, 
2202-2007 and 2358+0216; open circles in figure 4) have 
a peculiar continuum shape which neither of the empiri- 
cal extinction laws can explain. In this case astronomical 
silicate with small grain size of about 0.03 to 0.04 /im best 
reproduces the extinction in all three bands. 

The hi-lo and hi-BALs, on the other hand, have a 
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wider distribution near to the "origin" of the curves in the 
TCDD. However their small amount of extinction makes 
it difficult to uniquely distinguish the type of extinction. 

Table 1 (column 6) shows the determined extinction 
type for each object. For those objects which have posi- 
tions above the silicate-line (grain size 0.001 /zm) no ex- 
tinction type could be derived and they have been marked 
with "-" in table 1. 

Remarkable is that most of the observations in the 
TCDD do not lie on the curve produced with the MW 
extinction law. There are only two quasars that have 
extinction type LMC-MW in our sample. Since F2100 is 
very sensitive to the 2175 A feature, we can conclude from 
the lack of MW type extinction objects that the 2175 A 
bump is virtually absent in BAL QSOs. For the same 
reason, theoretical extinction curves of graphite grains 
(Laor, Drain 1993), which have pronounced features at 
around 2200 A, are unsuitable for dereddening of BAL 
QSO spectra. 

There is a clear division of the determined extinction 
laws (table 1) between lo-BALs and hi-BALs. This can 
be due to different viewing angles amongst various types 
of BAL QSOs (Yamamoto 1998) which means that the 
areas where the extinction occurs are in separate regions 
where different dust composition, grain size distribution, 
and physical conditions might affect the extinction law. 
It might also be due to a strong selection effect of the 
lo-BALs since only a few objects are known. 



4-.3. Extinction with Spatial Distribution of Dust 

The spatial distribution of dust is regarded as impor- 
tant for determining the effect of extinction on the spec- 
trum in the case of galaxies (e.g. Witt et al. 1992; Calzetti 
et al. 1994). In the case of quasars, however, there is 
only one continuum source with most of the continuum- 
forming region being too hot for dust to survive. There- 
fore, the uniform dust-screening model commonly used 
for stars might be the most suitable configuration of dust 
for investigations of extinction effects in quasars. 

To demonstrate the robustness of the TCDD in the 
case of more arbitrary spatial distribution of the dust in 
extended regions, we apply the attenuation curve (AC) 
from Calzetti et al. (1994) to FSED. Figure 4 illustrates 
that the AC-line lies on the left side close to the SMC-line 
and shows that a more complicated configuration than 
simple screening does not significantly affect the shape of 
the lines in the diagram. We expect that more realistic 
treatment of the dust distribution around quasars and 
solution of the radiative transfer problem are unlikely to 
endanger our assumption of the dust screen model where 
the absorption occurs outside the BAL region. 



4-. 4- Dereddening of Quasar Spectra 

The method employed by Sprayberry and Foltz (1992) 
to show evidence for extinction in BAL QSOs used only 
composite spectra from Weymann et al. (1991). We show 
in our diagram (figure 4) that the composite spectra they 
used, namely the hi- and lo-BALs, are both affected by 
dust. This means that estimation of the extinction in lo- 
BALs was done with respect to the intrinsically reddened 
composite spectrum of hi-BALs and can only be regarded 
as a hint for dust extinction in quasars. 

In our case the type and amount of intrinsic reddening 
of BAL QSOs are obtained from the TCDD by look- 
ing for the observed points' nearest positions to the lines 
produced by the extinction laws. Since we apply only 
standard extinction laws (SMC, LMC, MW, and astro- 
nomical silicate) we assign, for simplicity, one type for 
the cases where positions in the TCDD indicate an inter- 
mediate type of extinction. 

Several trials of the dereddening process were made 
for each object with the suggestion offered by the TCDD 
as first guess. In most cases, the extinction predicted 
by the TCDD gave satisfactory results. Then the spec- 
tra in the region around C IV and Si iv were examined 
for the quality of the fit and the amount of extinction 
slightly readjusted by eye to match the continuum level 
between C iv and Si iv. A more systematic fit would 
be possible if another band could be defined in this re- 
gion. However, the continuum in this region is subjected 
to the C IV broad absorption line, whose maximal width 
strongly varies from object to object. This makes it im- 
possible to define a band between Civ and Siiv that 
would serve for all BAL QSOs as a measure of the con- 
tinuum level. 

4-. 5. Examples of Dereddened Spectra 

Six objects are chosen (marked with asterisk in table 1) 
for presentation of the dereddened spectra as follows: 
the object with the highest extinction (1556-f3517) to 
show that the recovery of the spectrum works correctly 
for extreme cases; at intermediate extinction, representa- 
tives of various extinction types are chosen (2225—0534, 
0840+3633, and 1331—0108, which we have determined 
to be LMC, SMC, and silicate type respectively) to show 
that spectra dereddened with the predicted extinction 
types indeed recover the overall shape of the spectrum; 
two are taken from the region with the highest popula- 
tion density in TCDD (1228+1216 and 1303+3080) with 
one of the objects near to the MW-line in figure 4. 

Figure 5 shows dereddened and observed spectra to- 
gether with the corresponding power-law continua. Nu- 
merous emission lines significantly change the local shape 
of the quasar continuum which cannot be presented with 
the power-law and therefore we show FSED in panel (a) 
as a guide of the continuum level. 
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Fig. 5. Comparison of dereddened spectra (solid spectrum) with observed spectra (dotted line). The continua shown are the power-law 
representation of the Francis et al. (1991) composite spectrum without (upper smooth line) and with (lower smooth line) reddening. 
The spectra are normalized at 2100 A for the unreddened power-law. The object name, type of extinction law, and E{B — V) are 
indicated. The positions of the bands where the flux was measured are indicated as black rectangles at the bottom of each panel. 
Panels (a) and (b) show expamles for SMC type extinctions. In panel (a) of this figure we have over-plotted the composite spectrum 
from Francis et al. (1991). In panel (b) the observed spectrum and the reddened continuum have been scaled by a factor of 15 since 
they would be almost fiat in the present scale. 



The iron lo-BAL 0840+3633 in panel (a) was dered- 
dened with SMC type extinction. The successful recovery 
of the spectrum shows that strong Mg ii and Fe ii absorp- 
tion lines do marginally affect the measurements of the 
flux. It is striking to see that even the heavily extinct 
spectrum of 1556-1-3517 (panel b) can be recovered with 
the extinction values determined from the TCDD. The 
ratio of the flux in B2100 in this case between observed 
and dereddened spectrum is a factor of ~ 60. For these 
two objects adjustment of E(B — V) to values slightly 



smaller than suggested by the TCDD was necessary. 

The dereddened spectrum of the weakly extinct hi- 
BAL 1228+1216 (panel c) almost exactly matches FSED. 
The significantly stronger reddened hi-BAL 2225—0534 
(panel d) for which the LMC type extinction curve was 
applied also is in good agreement when compared with 
FSED. Note that FSED is only shown in panel (a) of 
figure 5. 

Dereddening of the hi-BAL 1303+3080 (panel e), 
which was classified as LMC-MW, with a MW type ex- 
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Fig. 5 (continued). Panels (c) and (d) show cxpamlcs for LMC type extinctions. 



tinction law gave somewhat better resuhs than did the 
LMC type. In the TCDD (figure 4) this object is closest 
to the MW-linc for which measurable extinction could be 
derived. We show one example for astronomical silicate 
type extinction from our sample — 1331—0108 (panel f). 
In this case, dereddening with the extinction values de- 
rived from the TCDD also leads to satisfactory results. 
We note that dereddening of 1331-0108 with SMC type 
extinction law, as figure 2 (panel a) would suggest, cannot 
recover the region around C iv whereas the steep change 
in the silicate extinction curve in this region gives a bet- 
ter fit. 

In general, better fits were achieved with the em- 
pirical extinction laws. The only exceptions from the 
whole quasar sample were 1331—0108, 2202—2007, and 



2358+0216, for which only astronomical silicate with a 
grain size smaller than for the other objects could re- 
cover the hump-like shape of the continuum at around 
1750 A. 

5. Conclusions 

26 BAL QSOs have been analyzed in this paper for 
which we have defined bands to measure the fiux and con- 
structed a two-colour diagnostic diagram (TCDD). Based 
on the TCDD we introduced a classification scheme for 
the determination of the extinction type and estimation 
of the amount of extinction for individual objects. Using 
this information we can deredden observed spectra with 
standard extinction laws (SMC, LMC, MW, and astro- 
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Fig. 5 (continued). Panel (e) shows an example for MW type extinction. In the case of astronomical silicate (f) the grain size in fim 
is given instead of the type of extinction law. 



nomical silicate). From our investigation we draw the 
following conclusions: 

For low-ionization BAL QSOs, the SMC type extinc- 
tion law is appropriate for dereddening while the range 
of extinction laws for high-ionization BAL QSOs varies 
from Milky Way, LMC to SMC. 

The TCDD also indicates that the Milky Way type 
extinction curve is not applicable for most of the BAL 
QSOs. 

The Francis ct al. (1991) composite spectrum is more 
suitable than Zheng et al. (1997) composite spectrum for 
the representation of the continuum shape of BAL QSOs. 

The similarity of the shapes of BAL QSOs' dered- 
dened continua to those of non-BAL QSOs suggests that 



both types of quasar have the same kind of continuum 
source. This means that a representative non-BAL (nor- 
mal) QSO spectrum such as the composite spectrum from 
Francis et al. (1991) might be appropriate for determining 
the continuum shape of incident radiation for the broad 
line region. 

There is a hint of the existence of an "intrinsic" con- 
tinuum shape of quasars that is bluer than the Francis 
et al. (1991) composite spectrum. This would imply that 
all quasars observed might have spectra that are red- 
dened by dust. 
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